Effect of a Normal-State Pseudogap on Optical Conductivity in Underdoped Cuprate 

Superconductors 
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We calculate the c-axis infrared conductivity a c (u>) in underdoped cuprate superconductors for 
spinfluctuation exchange scattering within the Cu02-planes including a phenomenological d-wave 
pseudogap of amplitude E g . For temperatures decreasing below a temperature T* ~ Eg/2, a gap 
for uj < 2E g develops in a c (u) in the incoherent (diffuse) transmission limit. The resistivity shows 
'semiconducting' behavior, i.e. it increases for low temperatures above the constant behavior for 
Eg = 0. We find that the pseudogap structure in the in-plane optical conductivity is about twice 
as big as in the interplane conductivity <7 c (w), in qualitative agreement with experiment. This is a 
consequence of the fact that the spinfluctuation exchange interaction is suppressed at low frequencies 
as a result of the opening of the pseudogap. While the c-axis conductivity in the underdoped regime 
is described best by incoherent transmission, in the overdoped regime coherent conductance gives a 
better description. 
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Numerous experiments have established the fact that 
the underdoped cuprate superconductors exhibit a 'pseu- 
dogap' behavior in both spin and charge degrees of free- 
dom below a characteristic temperature T* which can be 
well above the superconducting transition temperature 
T c . Many interpretations of the pseudogap have been 
advanced (see, e.g. the discussion in Ref. |l|), however, no 
consensus has been reached so far, which of the various 
microscopic theories is the correct one. It has been shown 
by Williams et aljj that specific heat, susceptibility and 
NMR data of many underdoped cuprates can successfully 
be modeled using a phenomenological normal-state pseu- 
dogap having <i-wave symmetry and an amplitude which 
is temperature independent but increases upon lowering 
the doping level into the underdoped regime. This strong 
anisotropy of the pseudogap is also in accordance with 
angle-resolved photoemission spectroscopy (ARPES) g 
periments on underdoped Bi 2 Sr 2 CaCu208-<5 (Bi2212)B' 
This model yields a smooth evolution of the normal state 
pseudogap into the superconducting gap as has been 
found in scanning-tunneling (STM) experimentsQ. Also, 
measurements of resistivity, Hall coefficient and,jthermo- 
electric power can be reconciled with this modellm. How- 
ever, measurements of the dynamical conductivity within 
or perpendicular to the Cu02-planes consistently show 
pseudogap structures having a size which only weakly 
depends on doping in the underdoped regime, in marked 
contrast to the model above. In addition, while the c-axis 
conductivity in YBa2Cu307_ d - (YBCO) shows a pgeiido- 
gap having a size of approximately 300-400 cm -1 .Em the 
size of the pseudogap extracted from ahjalane conduc- 
tivity is of the order of 600-700 cm^ffl This differ- 



ence cannot be attributed to the charge reservoir layers 
between the Cu02-planes, as it has been convincingly 
shown recently that the pseudogap seen in. c-axis con- 
ductivity has its origin in the CuCVplanesEJ. 

Here, we try to address these apparent inconsisten- 
cies in an effort to come to a consistent phenomenolog- 
ical description of the pseudogap. We calculate the c- 
axis and ab-plane conductivity in the presence of a tem- 
perature independent, but doping dependent pseudogap 
as is suggested by the work of Williams et a£|. To do 
so it is necessary to take into account some scattering 
mechanism within the Cu02-plane. We choose to study 
spinfluctuation exchange scattering within the so-called 
selfconsistent fluctuation-exchange (FLEX) approxima- 
tion for the two-dimensional Hubbard modellij, which 
has proven to give a good qualitative descriptiim-Qf the 
high-T c cuprates in the optimally doped regimeEj E£l. Es- 
pecially, the selfconsistently calculated interaction due 
to exchange of spin and charge fluctuations yields a 
quasiparticle scattering rate which varies linearly with 
frequency in the normal state, and exhibits a gap-like 
suppression at lower frequencies in the superconducting 
state. At the same time the effective mass ratio is in- 
creased at lower frequencies in the superconducting state. 
Thus, FLEX approximation accounts for the damping 
and mass enhancement needed to extend the theory for 
thermodynamical quantities by Williams et al.EJ to dy- 
namical quantities. 

We have shown previously that such a theory is capable 
of describing qualitatively a number of different experi- 
ments in underdoped cuprates like Knight-shift, nuclear- 
spin relaxation rate, neutron scattering, ARPES, STM, 
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and ab-plane conductivity experiment sl 16 li 17 l. In partic- 
ular, the pseudogap leads to a gap-like suppression of 
the scattering rate at lower frequencies below the linear 
frequency extrapolation in accordance with the observed 
in-plane conductivity. The T c which is calculated selfcon- 
sistcntly from the strong-coupling gap equation is sup- 
pressed in proportion to the magnitude of the pseudogap. 
The behavior of the Knight-shift, nuclear-spin relaxation 
rate, and density of states is described correctly as the 
temperature is decreased through T c , showing a smooth 
evolution of the pseudogap into the superconducting gap. 

Here, we will show that the selfconsistency of the spin- 
fluctuation interaction with the single particle proper- 
ties, especially with the pseudogap itself, provided by the 
FLEX approximation will lead to a natural understand- 
ing of the differences in size of the pseudogap seen in 
ab-plane and c-axis conductivity. Also, the semiconduct- 
ing behavior_af_the c-axis resistivity can be understood 
qualitativelyu'uEj. However, we will also see that such a 
simple ansatz for the pseudogap is not sufficient to un- 
derstand the doping independence of the size of the gap 
observed in ab-plane and c-axis conductivity. 

Measurements of the c-axis conductivity suggest that 
conductance in c-direction is coherent in the overdoped 
regimeliil successively becoming incoherent in the under- 
doped regimeBS. We will therefore study the two limits 
of coherent and incoherent c-axis conductivity. We will 
see that indeed within our model incoherent conductance 
gives a good description of the underdoped regime, while 
coherent conductance is more appropriate for the over- 
doped regime, confirming previous interpretations of the 
c-axis conductivity. 

The coherent conductivity along the interplane c- 
directioH, is given to lowest order in the inter-layer hop- 
ping t±!t3 by 
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where cq and a are the c-axis and ab-plane lattice con- 
stants. Here, TV is the normal spectral function and A\ 
and A g are the anomalous spectral functions with respect 
to the superconducting gap and the pseudogap, respec- 
tively. These spectral functions are taken from a selfcon- 
sistcnt solution of the FLEX equations in the presence of 
the pseudogap and are given by: 
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FIG. 1. The coherent dynamical c-axis conductivity (Eq. 
(0)) for three temperatures T=0.1t (solid line), 0.05t 
(dashed), and 0.03i (dotted) (a) for a pseudogap amplitude 
Eg = (b) for Eg = 0.15t The dashed-dotted line in (a) 
holds for the superconducting state for E g = (T c = 0.023t) 
at T = 0.0171 



We want to emphasize that it is necessary to include the 
bubble contribution due to A g into the conductivities and 
susceptibilities. Neglection of this term leads to severe 
disagreement with the data. Following Ref. |l|, for the 
pseudogap we assume the form 



Eg (k) = Eg [cos k x — cos k y ] 



(5) 



where E g is temperature independent and increases with 
decreasing doping level below the optimal doping level. 
The FLEX equations are solved selfconsistently in the 
presence of this pseudogap and yield the quasiparticle 
selfenergy components Z(k, w) and £(k, ui) as well as the 
superconducting gap 0(k, oj). Within FLEX the effective 
spin and charge fluctuation interactions are given by the 
RPA expressions 
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where the bubble spin susceptibility XsO is calculated self- 
consistently from the expression 
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ImXso(q, w) = 7T / dw' [/(a/) - f(u>' + ui)] 
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and thus depends on the pseudogap via the spectral func- 
tions. The charge susceptibility \c0 has the opposite sign 
in front of the anomalous terms in Eq. (Q) . For the on- 
site Coulomb repulsion we take an effective J7(q) with 
maximum value U = 3.6 at q = Q = (tt, 7r), as has been 
discussed in Ref. |l7]. For the bandstructure we take 
the tight-binding band 

ej; = t [—2 cos k x — 2 cos k y + AB cos k x cos k y — fi] (8) 

where t is the in-plane hopping matrix element and we 
take B = 0.45 and fi = —1.1 which describes approxi- 
mately the Fermi surfaces of the YBCO and Bi2212 com- 
pounds. 

In Fig. 0(a) we show the coherent conductivity <t c (uj) 
calculated in this way from Eq. (|l]) for different tempera- 
tures. Here we took E g = 0. For decreasing temperatures 
T a coherent Drude peak develops at low frequencies. 
Such a development of a coherent, Drude peak has been 
observed in overdoped cupratesOEl, where the pseudogap 
is absent or small. Thus, our coherent conductance re- 
sults account well for this observation in the overdoped 
compounds. In the superconducting state a suppression 
of a c {to) at intermediate frequencies sets in as shown by 
the dashed-dotted line in Fig. [l](a) for T = 0.017*. Here, 
T c = 0.023t. At the same time the Drude peak continues 
to sharpen. Fig. |l|(b) shows the normal-state coherent 
conductivity <J c {u>) in the presence of the pseudogap with 
amplitude E g — Q.lbt. While the pseudogap leads to a 
suppression of a c (u) at intermediate frequencies, the co- 
herent Drude peak at low frequencies still remains and 
even sharpens, similar to the superconducting state in 
Fig. |l](a). These results are completely different from 
the experimental results in the normal state of under- 
doped cuprates, showing instead of a coherent Drude 
peak a gap- like suppression at low frequencies. Thus, 
it is not sufficient to simply turn on a pseudogap in or- 
der to account for the c-axis conductivity, in underdoped 
compounds. As has been noted earlieraiiij the c-axis con- 
ductance at the same time becomes incoherent in the un- 
derdoped regime and therefore it is necessary to calculate 
the incoherent conductivity in the presence of a pseudo- 
gap. Incoherent conductivity corresponds to diffuse c- 
axis transmission and amounts in taking the averages of 
the spectral functions N(k,u)), Ai(k, w), and A g (k,uj) 
over all momenta (see the discussion in Ref. [D^). This 
means that N(k, uj) is replaced by the density of states 



symmetry of the superconducting and the pseudoga 
Then we find: 
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FIG. 2. The incoherent dynamical c-axis conductivity 
(Eq. @) for three temperatures T=0.1t (solid line), 0.05t 
(dashed), and 0.03t (dotted) and two values of E g =0 (upper 
three curves) and 0.15* (lower three curves). 

Fig. H shows the incoherent conductivity crj. ncoh (oj) for 
Eg = and E g = 0.15t for three different temperatures 
T = 0.1, 0.05, and 0.03*. For E g = 0.15* a gap de- 
velops below a threshold frequency of about uj ~ 2E g 
upon lowering the temperature while the conductivity 
stays almost constant for frequencies above this thres- 
hold energy. Here, E g = 2E g /KcZ(E g ) is the renormal- 
ized amplitude of the <i-wave pseudogap, and WeZ{ui) is 
the average mass renormalization at the Fermi surface 
which is of the order of 2 for the parameters considered 
here. For E g — 0, a" lcoh is very much frequency and 
temperature independent. These results are in qualita- 
tive agreement with the measured intgrplane conductiv- 
ity in underdoped YBCO compounds^. The gap in the 
c-axis conductivity a-J, ncoh (a;) for frequencies ui below 2E g 
develops below a characteristic temperature T* ~ E g /2. 
We want to stress that the temperature evolution of all 
physical quantities arises exclusively from the Fermi and 
Bose functions occuring iia— the FLEX equations in our 
real-frequency formulationO-3 and in the expressions for 
the susceptibilities (see Eq. @)), since we assumed that 
the pseudogap E g (]t) defined in Eq. (||) is temperature 
independent, following Ref. |l|. Physically, this means 
that above T* the effect of the pseudogap is smeared 
out such that the normal state behavior (corresponding 
to the FLEX equations for E g = 0) is recovered while 
the effect of the pseudogap on the quasiparticle and spin 
excitation spectra increases as T decreases below T* to- 
wards T r . 

From our incoherent conductivity we can extract the 
while the averages of A\ and A g vanish due to the d- wave c-axis resistivity p c — [a" icoh (uj — 0)]" 1 in the presence 
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FIG. 3. The incoherent c-axis resistivity p c as a function of 
temperature for E g =0 (solid line), O.lt (dashed), 0.15t (dot- 
ted), and 0.2t (dashed-dotted). 



of the pseudogap. The temperature dependence of p c is 
shown in Fig. || for different values of the pseudogap 
amplitude E g =0, 0.1, 0.15, and 0.2t, corresponding to 
different doping levels. For E g = 0, p c is almost constant. 
For finite E g it starts to increase above the curve for E g 
at lower temperatures. This 'semiconducting' behavior 
of p c is directly related to the depth of the pseudogap 
at zero frequency in a c ncoh (ui = 0) which increases for 
increasing E g and decreasing temperature (see Fig. ||). 
The results in Fig. ^ are consistent with the estimate 
of the characteristic temperature given above, i.e. T* ~ 
E g /2 = E g /ReZ, because the steep rise of p c for a given 
E g appears approximately below T*. Our results are 
in qualitative agreenient with c-axis resistivity data on 
underdoped YBCOola. We remark that the definition of 
the characteristic temperature T* ~ E g /2 corresponds to 
the scaling procedure in Ref. where it has been shown 
that the NMR Knight shift ™K n (T) for a wide range 
of doping values follows closely a universal scaling curve 
if the data are plotted against a scaling parameter z = 
2T/E g . The downturn of S9 K n (T) for decreasing z occurs 
at about z = 1 which corresponds to a T* ~ E g /2. 

In Fig. |i| we show a c ncoh (uj) at a fixed temperature 
T = 0.03t and different values of E g . From Fig. ^ we see 
that the renormalized size 2E g of the gap in crJ. ncoh follows 
E g . If one assumes that E g strongly changes with doping 
level, as has been proposed in Ref. [l], the results shown in 
Fig. U cannot provide an explanation for the doping in- 
dependence of the pseudogap seen in c-axis conductivity 
on underdoped cuprates. This is an apparent inconsis- 
tency of the model by Williams et al, which can describe 
well thermodynamic quantities, but does not give a sat- 
isfactory account of the doping dependence of ab-plane 
and c-axis conductivity, which are dynamical quantities. 
Our results rather indicate that there are two indepen- 
dent energy scales involved in the pseudogap problem: 
first, the width of the pseudogap (here 2E g ) as it is seen 
in the w-dependence of the c-axis and ab-plane conduc- 
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FIG. 4. The incoherent dynamical c-axis conductivity (Eq. 
@) for T = 0.03t and E g =G (solid line), O.lt (dashed), 
0.15i (dotted), and 0.2i (dashed-dotted). For comparison we 
also show the result in the superconducting state for E g — 
(T c = 0.0234) at T — 0.017* (lower solid line). 



tivity, being largely doping independent, and second, the 
depth of the pseudogap as it is seen in the u> = value 
of the incoherent c-axis conductivity corresponding to 
the characteristic temperature T* for the c-axis resistiv- 
ity and the thermodynamic quantities, which increases 
upon lowering the doping level. Such two energy scales 
could be introduced into the problem by considering more 
complicated forms for the pseudogap than Eq. (||). For 
example, the pseudogap could have a frequency depen- 
dence or a momentum dependence which changes with 
temperature as suggested by a recent analysis of ARPES 
datJM 

The lower solid line in Fig. [| shows the result for 
the incoherent c-axis conductivity in the superconduct- 
ing state for T = 0.017* and E g = (T c = 0.023t). This 
compares wellrwith the experimental results on optimally 
doped YBCOjD showing a suppression at low frequencies, 
similar to the suppression due to the normal-state pseu- 
dogap. In addition, a weak enhancement at the gap edge 
develops. Again, this behavior is completely different 
from the corresponding behavior of the coherent c-axis 
conductivity (see the dashed-dotted line in Fig. |l|(a) ). 

The in-plane conductivity a a b(uj), neglecting vertex 
corrections, is given by 

<7^ («) = ■= / d/[/(a/)-/(u/ + W )] 

k 

+ Ai(k,cj' + w)Ai(kV) + A g (k,uj' + w)A g (k,u')] , 

where Vk,i = dek/dki are the band velocities within the 
ab-plane. Our results for a a b(to) (see Fig. 6 in Ref. [l7]) 
are very similar to those for the coherent c-axis conduc- 
tivity shown in Fig. |l|. The in-plane conductivity is co- 
herent in character and shows a Drude peak at low fre- 
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quencies even in the underdoped compounda3'c2l. How- 
ever, the size of the pseudogap structure seen in ab-plane 
conductivity on underdoped YBCO has been measured 
to be 600-700 cm -1 ,BEJ while the gap extracted from c- 
axis conductivity in. the same compounds is of the order 
of 300-400 cm' - 1 J3'H In Fig. | we show our results for 
CTah(w) and a™ coh ((j) for E g = 0.15* and T = 0.03< along 
with the density of states N(oj) Eq. (^|) and the quasi- 
particle damping rate admZ(k a , uj) at the antinodal mo- 
mentum k a at the Fermi surface (arbitrary units). Here 
we see that the size of the pseudogap appearing in these 
four quantities is quite different. In fact, the gaps have an 
approximate relation of 1:2:3:4 in the density of states, in- 
coherent c-axis conductivity, quasiparticle damping rate, 
and ab-plane conductivity, respectively. Especially, we 
find that the gap structure in <j a b(u)) is about twice as 
big as in crJ, ncoh (w), being about 4E g in a a b(oj) while only 
2E g in cr" lcoh (w) with E g ~ 0.12i, in rough agreement 
with experiment. This relation of the gaps is a direct 
consequence of the electronic origin of the spinfluctuation 
scattering process and the selfconsistency of the FLEX 
equations: the opening of the pseudogap leads to a sup- 
pression of the spinfluctuation interaction via Eqs. (Q) 
and (^) at frequencies below ~ 2E g . This in turn results 
in a reduction of the sclf-cncrgy below ~ 3E g and a corre- 
sponding structure at ~ AE g in a a ),(u)). However, in the 
incoherent c-axis conductivity only a gap of size ~ 2E g 
appears because of the momentum average of the spec- 
tral functions, resulting in a frequency convolution of the 
density of states with itself. The appearance of a 4Ao- 
gap in the ab-plane conductivity in the superconduct- 
ing state for an electronic pairing mechanism has been 
noted .eaxlier in connection with marginal Fermi liquid 
theorycirEl Here we suggest that a corresponding ef- 
fect is taking place in the pseudogap state of underdoped 
high-T c compounds. 

The pseudogap structures in the curves in Fig. [5] ap- 
pear to be washed out somewhat and show more com- 
plex behavior than a simple suppression at the effective 
pseudogap. This is due to the fact that the pseudogap 
Eq. (||) is renormalized due to self-energy effects. The 
structures seen in the conductivity, density of states, and 
quasiparticle damping rate do not display a pure <i-wave 
gap, but a renormalized one, similar as in the supercon- 
ducting state (see Rcf. [is]). More detailed discussions of 
the density of states and the quasiparticle damping rate 
can be found in Refs. [if] [IT]. 

To summarize, we have investigated the influence 
of a normal-state pseudogap of the form suggested by 
Williams ct allil on the c-axis and ab-plane conductivity 
for spinfluctuation exchange scattering within the self- 
consistent FLEX approximation. We find that coher- 
ent conductance can describe the c-axis conductivity in 
the overdoped compounds, while it is necessary to con- 
sider incoherent c-axis conductance in the underdoped 
regime. Incoherent conductance can account well for the 
dynamical c-axis conductivity and the c-axis resistivity 
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FIG. 5. Density of states N(ui) (solid line), incoherent 
c-axis conductivity <7j, ncoh (a;) (dashed), quasiparticle damping 
rate wImZ(k a ,[j) (dotted), and ab-plane conductivity cr a b(ui) 
(dashed-dotted) as a function of frequency for E g — 0.15i 
and T = 0.03t (arbitrary units). All four quantities show 
gap-like suppressions at low frequencies. The sizes of these 
gaps roughly have a relation of 1:2:3:4. The ab-plane conduc- 
tivity a a b{io) shows a strong Drude peak at low frequencies 
within the gap. 

in the underdoped compounds, showing 'semiconducting' 
behavior. However, it is difficult to reconcile the doping 
dependence of the amplitude of the pseudogap, as sug- 
gested by the work of Williams et al, with the doping 
independent size of the pseudogap seen in dynamical c- 
axis and ab-plane conductivity. This suggests that the 
pseudogap has a nontrivial momentum or frequency de- 
pendence, which changes with temperature. We find that 
the difference in size of the pseudogap in ab-plane con- 
ductivity as opposed to c-axis conductivity finds a natu- 
ral explanation in the electronic origin of spinfluctuation 
scattering and its selfconsistency with the single-particle 
properties. This leads to a gap structure of size ~ 4E g 
in the ab-plane conductivity, while the gap seen in the 
incoherent c-axis conductivity only has a size of ~ 2E g . 
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